Dear all,
I have a matrix of reads counts for each gene for each strain of an animal. I am using DESeq2 to analyse the clustering of the strains. I have performed a Variance Stabilising Transform on my DESeq2 object and...can do a PCA:
p <- DESeq2::plotPCA(vsd, intgroup=c('Condition'))
p <- p + labs(title = 'PCA on gene expression levels')
p
Now, I would like to know the coeffi…